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Various recent studies proved that cosmological models with a significant contribution from cold 
dark matter isocurvature perturbations are still compatible with most recent data on cosmic mi- 
crowave background anisotropies and on the shape of the galaxy power spectrum, provided that one 
allows for a very blue spectrum of primordial entropy fluctuations (ni so > 2). However, such models 
predict an excess of matter fluctuations on small scales, typically below 40/i -1 Mpc. We show that 
the proper inclusion of high-resolution high signal-to-noise Lyman-a forest data excludes most of 
these models. The upper bound on the isocurvature fraction a = f? so / (1 + /j| ), defined at the pivot 
scale fco = 0.05 Mpc -1 , is pushed down to a < 0.4, while rii so = 1.9±1.0 (95% confidence limits). We 
also study the bounds on curvaton models characterized by maximal correlation between curvature 
and isocurvature modes, and a unique spectral tilt for both. We find that /i so < 0.05 (95% c.l.) in 
that case. For double inflation models with two massive inflatons coupled only gravitationally, the 
mass ratio should obey 7? < 3 (95% c.l.). 
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I. INTRODUCTION 

With the most recent measurements of the cosmic mi- 
crowave background (CMB) anisotropies and large scale 
structures (LSS) of the universe as well as various other 
astronomical observations, it is now possible to have a 
clear and consistent picture of the history and content of 
the universe since nucleosynthesis. In particular, it is well 
established that the cosmological perturbations which 
gave rise to the CMB anisotropies and the LSS of the 
universe were inflationary, with a close to scale-invariant 
Harrison-Zeldovich spectrum. Moreover, the CMB and 
LSS data allow to test the paradigm of adiabaticity of the 
cosmological perturbations and hence the precise nature 
of the mechanism which has generated them. 

The simplest realizations of the inflationary paradigm 
predict an approximately scale invariant spectrum of 
adiabatic (AD) and Gaussian curvature fluctuations, 
whose amplitude remains constant outside the horizon, 
and therefore allows cosmologists to probe directly the 
physics of inflation from current CMB and LSS ob- 
servations. However, this is not the only possibility. 
Models of inflation with more than one field generi- 
cally predict that, together with the adiabatic compo- 
nent, there should also be entropy, or isocurvature per- 
turbations P, 0, Oli EH- [3- E3l - associated with fluctuations 
in number density between different components of the 
plasma before photon decoupling, with a possible statis- 
tical correlation between the adiabatic and isocurvature 
modes 0. Baryon isocurvature (BI) perturbations and 
cold dark matter isocurvature (CDI) perturbations were 
proposed long ago Q as an alternative to adiabatic per- 
turbations. These BI and CDI modes are qualitatively 
similar, since they are related by a simple rescaling factor 
Qg/f^ dm , or ilB/^cdm for the cross-correlation: thus, by 



studying the case of mixed AD + CDI modes, one implic- 
itly includes the case of AD + BI, for which the allowed 
isocurvature fraction is larger roughly by the above fac- 
tor evaluated near the maximum likelihood model. A few 
years ago, two other modes, neutrino isocurvature den- 
sity (NID) and velocity (NIV) perturbations, have been 
added to the list Moreover, isocurvature perturba- 
tions have been advocated in order to explain the high 
redshift of reionization claimed by the WMAP team [l(| ■ 

Note, however, that in the case all fields thermalize at 
reheating, no isocurvature mode will survive [Til ]. The 
simplest assumption for generating observable CDI per- 
turbations is that one of the infiaton fields remains un- 
coupled from the rest of the plasma between inflation 
and its decay into CDM particles. Since baryons and 
neutrinos arc usually assumed to be in thermal equilib- 
rium in the early Universe, it is more difficult to build 
realistic models for the generation of BI, NID and NIV 
modes than for CDI - but some possibilities still exist, 
based on non-zero conserved quantities and chemical po- 
tentials (see e.g. 

Moreover, it is well known that entropy perturbations 
seed curvature perturbations outside the horizon J-jUjf 
so that it is possible that a significant component of 
the observed adiabatic mode could be maximally cor- 
related with an isocurvature mode. Such models are 
generically called curvaton models 0, HdL [Til ITt3 | , and 
are now widely studied as an alternative to the standard 
paradigm. Furthermore, isocurvature modes typically in- 
duce non- Gaussian signatures in the spectrum of primor- 
dial perturbations |l7j . 

In the last few years, various models with a corre- 
lated mixture of adiabatic and isocurvature perturba- 
tions have been tested by several authors, with differ- 
ent combinations of data sets and theoretical priors. A 
crucial difference between these analyses lies in the as- 
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sumptions concerning the scale-dependence of the var- 
ious modes. Some groups assumed for simplicity that 
the adiabatic and isocurvature mode shared exactly the 
same scale-dependence 0, ^| , but enriched the anal- 
ysis by considering the full mixtures of several modes 
at a time (AD, CDI, NID, NIV). Other groups concen- 
trated on the (correlated) mixture of two modes only 
(AD+CDI in [21 EHH, AD+NID and AD+NIV in 
|21j), with a different power law for the three compo- 
nents (adiabatic, isocurvature and cross-correlation), as 
expected in the general case. Finally, an intermediate 
approach consists in studying the mixture of two modes 
with a scale- independent mixing angle, i.e., only two tilts 
H HI 0- In addition to these references, 
some groups studied the case of the curvaton scenario, 
which requires some specific analyses |2^, |2^, |2^| 
since it involves a maximal correlation/anti-correlation 
and a unique spectral index for the adiabatic and isocur- 
vature modes. Furthermore, two groups have quantified 
the need for isocurvature modes through a Bayesian Ev- 
idence computation on the basis of current CMB and 
galaxy power spectrum data, reaching somewhat differ- 
ent conclusions due to a different choice of priors (3^, El • 

In this work, we are particularly interested in mixed 
models with AD+CDI modes and three different tilts, 
for which it was shown in Refs. (2j| and that a 
significant fraction of isocurvature perturbations is still 
allowed. This sounds surprising at first sight, since 
the isocurvature mode is known for suppressing small- 
scale CMB anisotropics. This is true indeed for a scale- 
invariant spectrum of primordial isocurvature fluctua- 
tions, but not in general: a significant isocurvature con- 
tribution with a very blue tilt (nj SO ~ 3) can contribute to 
CMB anistropies even on small scales, and can be com- 
patible to some extent with the CMB temperature and 
temperature-polarization spectra, in spite of the small 
shift induced in the scale of the acoustic peaks. These 
models predict gencrically an excess of matter fluctua- 
tions on small scales. Using the shape and amplitude of 
the linear power spectrum derived from galaxy surveys 
at wavenumbers k < 0.15ft-/Mpc, one can exclude such 
an excess for wavelengths A = 2ir/k larger than 40 ft- -1 
Mpc. The main goal of this work is to push the con- 
straints down by making use of Lyman-a forest data, 
which probe large-scale structure at redshift z ~ (2 — 3) 
and on scales (1 — 40) /i -1 Mpc, in the mildly non-linear 
regime. Therefore, in any comparison between Lyman-a 
observations and linear theoretical predictions, it is nec- 
essary to take into account the non-linear evolution with 
N-body or hydrodynamical simulations. 

Usually, these simulations are carried under the as- 
sumption of adiabaticity. However, it is not difficult 
to generalize them to the case of mixed adiabatic plus 
isocurvature models. During matter domination, the 
perturbations seeded by each of the two modes are in- 
distinguishable: the only difference lies in their scale- 
dependence, but not in their nature or time-evolution. 
So, a given mixed model is entirely specified by a sin- 



gle matter transfer function, defined for instance soon 
after the time of equality. Therefore, the Lyman-a for- 
est data can be safely applied to non-adiabatic models 
provided that one takes into account the fact that the 
matter transfer function has more freedom than in the 
purely adiabatic case. In the following analysis, we will 
carefully take this point into account. 

We will use here the linear matter power spectrum in- 
ferred from two large samples of quasar (QSO) absorption 
spectra |32l J33j using state-of-the-art hydrodynamical 
simulations [34( combined with cosmic microwave back- 
ground data from the WMAP satellite [35J; as well as 
from the small-scale temperature anisotropy probed by 
VSA H3, CBI [13 and ACBAR [H; from the matter 
power spectrum measured by the 2-dcgrcc-Field Galaxy 
Redshift Survey (2dFGRS) and the Sloan Digital Sky 
Survey (SDSS) [4(|; and finally from the recent type la 
Supernova (SN) compilation of Ref. ■ We note that 
the cosmological parameters recovered from the data sets 
used in this paper are in good agreement with subsequent 
studies made by the SDSS collaboration using a differ- 
ent data set and a very different theoretical modelling 
(ElEHEEni])- This demonstrates that the analy- 
sis of the Lyman-a forest QSO spectra is robust and that 
many systematic uncertainties involved in the measure- 
ment are now better understood than a few years ago. 

The plan of the paper is as follows. In section II we de- 
scribe the notations we used for the isocurvature sector. 
In section III we introduce the Lyman-a data that we are 
employing. In section IV we discuss the general bounds 
on our full AD+CDI parameter space from Lyman-a, 
CMB, LSS and SN data using a Bayesian likelihood anal- 
ysis. We also check explicitly with a hydrodynamical 
simulation the robustness of our Lyman-a data-fitting 
procedure, and we address the subtle issue of the role of 
paramctrizations and priors on the isocurvature bounds 
and in the interpretations of our results. We also discuss 
the specific curvaton models with maximal anticorrela- 
tion and equal tilts for both adiabatic and isocurvature 
modes, as well as bounds on double inflation models. In 
section V we draw our conclusions. 



II. MIXED ADIABATIC/ISOCURVATURE 
MODELS 

A. Primordial spectra 

For the theoretical analysis, we will use the notation 
and some of the approximations of Ref. plj . During 
inflation, more than one scalar field could evolve suffi- 
ciently slowly that their quantum fluctuations perturbed 
the metric on scales larger than the Hubble scale during 
inflation. These perturbations will later give rise to one 
adiabatic mode and several isocurvature modes. We will 
restrict ourselves here to the situation where there are 
only two fields, <pi and 4>2, and thus only one isocurvature 
and one adiabatic mode. Introducing more fields would 
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complicate the inflationary model and even then, it would 
be rather unlikely that more than one isocurvature mode 
contributes to the observed cosmological perturbations. 

Therefore, the two-point correlation function or power 
spectra of both adiabatic and isocurvature perturbations, 
as well as their cross-correlation, can be parametrized 
with three power laws, i.e. three amplitudes and three 
spectral indices, 
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Here, 1Z stands for the curvature perturbation, and <S = 
(4dm — 3<5 7 /4) for the CDI perturbation. Both are eval- 
uated during radiation domination and on super-Hubble 
scales. We also introduced an arbitrary pivot scale ko, 
at which the amplitude parameters are defined through 
A = (^ ad ) 1/2 and B = [S 2 ^) 1 / 2 . In addition to the fact 
that curvature and entropy perturbations are generally 
correlated at the end of inflation, some extra correla- 
tion can be generated later by the partial conversion of 
isocurvature into adiabatic perturbations. The correla- 
tion angle A(/c) is in general a function of k, and in the 
above definitions, we approximated cos A(/c) by a power 
law with amplitude cos Afe and tilt n cor . So, we assumed 
implicitly that the inequality 



k 

IcosAfcJ ( — 
k 



< 1 



(2) 



holds over all relevant scales. We will enforce this condi- 
tion in the following analysis. 



B. CMB anisotropy power spectra 

The angular power spectrum of temperature and 
polarization anisotropics seen in the CMB today can 
be obtained from the radiation transfer functions 
for adiabatic and isocurvature perturbations, Of d (fc) 
and OJ so (fc), computed from the initial conditions 
(^rad(fc) Aad(fc)) = (1,0) and (0,1), respectively, and 
convolved with the initial power spectra, 
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Then, the total angular power spectrum reads 



Ci = a' c, 
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In many works (see for instance |23J,|2J|), the following 
parametrization is employed: 

d = A 2 [C'f + f? so Cf° + 2/ iso cos A k0 Cf 01 } , (4) 

where /j so = B/A represents the entropy to curvature 
perturbation ratio during the radiation era at k = k^. 
We will use here a slightly different notation, used before 
by other groups 



a sligntly 
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Ci = (A 2 +B 2 ) (1 - a) Cf d + a Cf ° + 2/3 



a)Cr 
(5) 

where a = B 2 /(A 2 + B 2 ) represents the isocurvature 
fraction at ko, and runs from purely adiabatic (a = 0) to 
purely isocurvature (a = 1), while (3 defines the correla- 
tion coefficient at ko, with (3 = +1(— 1) corresponding to 
maximally correlated(anticorrelatcd) modes. There is an 
obvious relation between both parametrizations: 



« = /iL/a + fL) , 



j3 = cos A^ 



(G) 



This notation has the advantage that the full parameter 
space of (a, 2/?y / a(l — a)) is contained within an ellipse. 
The North and South rims correspond to fully correlated 
(/3 = +1) and fully anticorrelated ((3 = —1) perturba- 
tions, with the equator corresponding to uncorrelated 
perturbations ((3 = 0). The East and West correspond to 
purely isocurvature and purely adiabatic perturbations, 
respectively. Any other point within the ellipse is an ar- 
bitrary admixture of adiabatic and isocurvature modes. 

We should emphasize that the three amplitude param- 
eters (A 2 + B 2 ), a and (3 are defined at k = ko, and that 
comparing bounds from various papers is straightforward 
only when the pivot scale is the same. For instance, in 
the simple case where n ac i = n; so , a is independent of kg, 
but this is not the case for (3: if n cm > 0, points within 
the (a, 2(3yja{\ — a)) ellipse are shifted vertically to- 
ward the edges of the ellipse when one increases ko and 
shifted toward the horizontal (3 = line when one de- 
creases fco. When n ac i ^ n- iso , both a and (3 depend on 
the pivot scale. In addition, by changing the prior on the 
amplitudes, a shift in the pivot scale affects the n; so like- 
lihood quite dramatically |22j . Throughout this paper, 
we will use ko = 0.05 Mpc -1 , which is the most frequent 
choice in the literature. This value corresponds roughly 
to the multipolc number Iq ~ 400. Therefore, the ra- 
tio C^o/Cfoo is roughly independent of the tilt values. 
For cosmological parameters close to the best-fit ACDM 
model, one finds Q00/C400 ~ 0.01. The smallness of 
this number comes from the fact that Q\ so (k) is strongly 
suppressed with respect to 0^ d (fc) for large wavenumbcrs. 
Indeed, the metric perturbations induced by isocurvature 
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perturbations remain small during radiation domination: 
so, for small scales entering early inside the Hubble ra- 
dius, the amplitude of the photon acoustic oscillations is 
also small (as can be seen via its transfer function). As 
a consequence, even if during radiation domination one 
has «S ra d(fco) ~ T^rad(fco) (corresponding to / iso ~ 1 or 
a ~ 0.5) the isocurvature mode contributes only to 1% 
of the observed anisotropy near Iq. Of course, if n; so is 
very different from n a d, there could still be a large isocur- 
vature contribution at either larger or smaller scales. 



scales (typically for wavenumbcrs k > 0.1/i/Mpc). We 
illustrate this behavior on Fig. for a particular set of 
AD+CDI models with two different values of the isocur- 
vature tilt, ?ii sci = 3 or rii SO = 2.2, and many possible 
values of (a, (3). The impact of the non-adiabatic contri- 
bution consists either in a smooth change of the effective 
slope on small scales, or more radically in a sharp feature 
(a pronounced break or a dip). The second situation can 
occur on relevant scales for large positive /3, and of course 
large enough values of a and rii so . 



C. Matter power spectrum 

Since in the following we will focus on the constraints 
induced on mixed AD+CDI models by the Lyman-a 
data, let us give a few details on the shape of the lin- 
ear matter power spectrum 

P(k) = (A 2 + B 2 ) [(1 - a) P ad (fc) + a P iso (k) 

+2P^a(l-a)P cor (k)}. (7) 

Here P ad and P lso arc computed from the initial condi- 
tions (72. r ad(fc),<S rad (/c)) = (1,0) and (0,1) respectively, 
exactly like 9 ad (fc) and <df°(k), and the cross-correlated 
term is simply given by 



P COT (fc) = -(fc/fc )™-'[P ad (fc)P IS °(fc)]5 



(8) 



where the minus sign comes from the fact that with our 
definition of S, a positive correlation (72. ra d<Srad) > 
in the early Universe implies a reduction of the matter 
power spectrum today, and vice-versa. 

In the limit k 3> k eq , where k eq corresponds to modes 
crossing the Hubble length at the time of equality, it is 
well known (see e.g. ^(|) that the power spectra obey, to 
first approximation, 

P ad (fc) oc (k/k y^-Hn(k/k eq ) 2 , (9) 
P iso (fc) cx (k/ko) n '°°-\ (10) 

which shows that for ni so ~ n a d the isocurvature con- 
tribution to the small-scale power spectrum is generi- 
cally much redder than the adiabatic one. The relative 
amplitude depends on the cosmological parameters. In 
the vicinity of the concordance ACDM model, one finds 
P iso (fc )/P ad (fc ) - 4 x 10" 3 for CDI. So, like for CMB 
anisotropics, we see that even when <S ra d(fco) ~ 7£ ra d(fco) 
in the early universe (i.e. /i so ~ 1 or a ~ 0.5), the isocur- 
vature contribution to the currently observed power spec- 
trum is only of the per cent order, at least near the pivot 
scale. However, for large rii so , the contribution may be 
large on small scales. 

Indeed, a large portion of the parameter region allowed 
by previous studies corresponds to a significant isocurva- 
ture fraction a > 0.1 and to a very blue tilt n; so > 1.5. In 
this case, the matter power spectrum is affected or even 
dominated by the non-adiabatic contribution on small 



III. PROBING THE MATTER POWER 
SPECTRUM WITH THE LYMAN-a FOREST IN 
QSO ABSORPTION SPECTRA 

It is well established by analytical calculation and hy- 
drodynamical simulations that the Lyman-a forest blue- 
ward of the Lyman-a emission line in QSO spectra is 
produced by the inhomogeneous distribution of a warm 
(~ 10 4 K) and photoionized intcrgalactic medium (IGM) 
along the line of sight. The opacity fluctuations in the 
spectra arise from fluctuations in the matter density and 
trace the gravitational clustering of the matter distribu- 
tion in the quasi- linear regime |47j . The Lyman-a forest 
has thus been used extensively as a probe of the matter 
power s pec trum on comoving scales of (1 — 40) h~ 1 Mpc 

111 IM III III 

The Lyman-a optical depth in velocity space u (km/s) 
is related to the neutral hydrogen distribution in real 
space as (see e.g. Ref. H3) : 



t(u) 



dyn m {y) V[u-y- v\\(y), b(y)] dy , 

(11) 

where ero,a = 4.45 x 10 18 cm 2 is the hydrogen Lya cross- 
section, y is the real-space coordinate (in km s" 1 ), V 
is the standard Voigt profile normalized in real-space, 
b = (2fcsT/mc 2 ) 1 / 2 is the velocity dispersion in units of 
c, H(z) the Hubble parameter, njji is the local density 
of neutral hydrogen and vu is the peculiar velocity along 
the line-of-sight. The density of neutral hydrogen can 
be obtained by sol ving the photoionization equilibrium 
equation (see e.g. |50f)- The neutral hydrogen in the 
IGM responsible for the Lyman-a forest absorptions is 
highly ionized due to the metagalactic ultraviolet (UV) 
background radiation produced by stars and QSOs at 
high rcdshift. This optically thin gas in photoionization 
equilibrium produces a Lyman-a optical depth of order 
unity. 

The balance between the photoionization heating by 
the UV background and adiabatic cooling by the expan- 
sion of the universe drives most of the gas with 5b < 10, 
which dominates the Lyman-a opacity, onto a power-law 
density relation T = To (1 + J;,) 7-1 , where the param- 
eters To and 7 depend on the reionization history and 
spectral shape of the UV background and 8b is the local 
gas overdensity (1 + Sb = pb/pb)- 
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FIG. 1: (Left) Matter power spectrum for a family of mixed AD+CDI models, with all parameter fixed except a and ft (the 
global normalization also varies in order to mantain a fixed amplitude on large scales). In particular, in all models we kept 
"-ad = 0.95, ni so = 3 and n cor = 0. The thick line stands for the pure adiabatic case (a — 0). The thin solid (red) lines show 
uncorrelated models (/3 = 0) with a = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. The lower blue (upper green) dashed lines show the maximally 
correlated models with j3 = 1 (anti-correlated with (5 = —1) for the same values of a. (Right) Same as the left plot, but with 
7ii so reduced to 2.2. 



The relevant physical processes can be readily mod- 
elled in hydrodynamical simulations. The physics of a 
photoionized IGM that traces the dark matter distri- 
bution is, however, sufficiently simple that considerable 
insight can be gained from analytical modeling of the 
IGM opacity based on the so called Fluctuating Gunn 
Peterson Approximation neglecting the effect of peculiar 
velocities and the thermal broadening [HI]]. The Fluc- 
tuating Gunn Peterson Approximation makes use of the 
power-law temperature density relation and describes the 
relation between Lyman-a opacity and gas density (see 
|33l.l52l ]) along a given line of sight as follows, 

Hz) <x (1 + S^)f T-<"(2) = A(z) (1 + h(z)T , (12) 

y ' v 3.5 J v 0.02 J Vsoook; 
(_h_y l [ H(z)/H y 1 ( r ffl y 1 

\0.65J V 3 - 68 / V 1 - 5 x 10-12 s V ' 

where /3 = 2-0.7 (7- 1) in the range 1.6-1.8, r ffl the HI 
photoionization rate, Hq = h 100 km/s/Mpc the Hubble 
parameter at redshift zero. For a quantitative analysis, 
however, full hydrodynamical simulations, which prop- 
erly simulate the non-linear evolution of the IGM and its 
thermal state, are needed. 

Equations l|ll[) and (|12f> show how the observed flux 
F = cxp (— t) depends on the underlying local gas den- 
sity pb, which in turn is simply related to the dark matter 
density, at least at lar ge s cales where the baryonic pres- 
sure can be neglected j53|. Statistical properties of the 
flux distribution, such as the flux power spectrum, arc 
thus closely related to the statistical properties of the 
underlying matter density field. 



A. The data: from the quasar spectra to the flux 
power spectrum 



The power spectrum of the observed flux in high- 
resolution Lyman-a forest data provides meaningful con- 
straints on the dark matter power spectrum on scales 
of 0.003 s/km < k < 0.03 s/km, roughly corresponding 
to scales of (1 — 40) ft, _1 Mpc (somewhat dependent on 
the cosmological model). At larger scales the errors due 
to uncertainties in fitting a continuum (i.e. in removing 
the long wavelength dependence of the spectrum emit- 
ted by each QSO) become very large while at smaller 
scales the contribution of metal absorption systems be- 
comes dominant (see e.g. (3^, Hil). In this paper, we will 
use the dark matter power spectrum that Viel, Haehnelt 
& Springel [34| (VHS) inferred from the flux power spec- 
tra of the Croft et al. [H (C02) sample and the LUQAS 
sample of high-resolution Lyman-a forest data [5^ • The 
C02 sample consists of 30 Keck high resolution HIRES 
spectra and 23 Keck low resolution LRIS spectra and 
has a median redshift of z = 2.72. The LUQAS sample 
contains 27 spectra taken with the UVES spectrograph 
and has a median redshift of z = 2.125. The resolution 
of the spectra is 6 km/s, 8 km/s and 130 km/s for the 
UVES, HIRES and LRIS spectra, respectively. The S/N 
per resolution element is typically 30-50. Damped and 
sub-damped Lyman-a systems have been removed from 
the LUQAS sample and their im pac t on the flux power 
spectrum has been quantified by [33j. Estimates for the 
errors introduced by continuum fitting, the presence of 
metal lines in the forest region and strong absorptions 
systems have also been made [3^, H3, 13 El • 
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B. From the flux power spectrum to the linear 
matter power spectrum 

VHS have used numerical simulation to calibrate the 
relation between flux power spectrum and linear dark 
matter power spectrum with a method proposed by C02 
and improved by and VHS. A set of hydrodynamical 
simulations for a coarse grid of the relevant parameters 
is used to find a model that provides a reasonable but 
not exact fit to the observed flux power spectrum. Then, 
it is assumed that the differences between the model and 
the observed linear power spectrum depend linearly on 
the matter power spectrum. 

The hydrodynamical simulations are used to determine 
a bias function between flux and matter power spectrum: 
Pp(fc) = b 2 F (k) P(k), on the range of scales of interest. 
In this way the linear matter power spectrum can be 
recovered with reasonable computational resources. |6fij 
This method has been found to be robust provided the 
systematic uncertainties arc properly taken into account 
|34l 1571 ] . Running hydrodynamical simulations for a fine 
grid of all the relevant parameters is unfortunately com- 
putationally prohibitive (see discussion in [4^ | on a pos- 
sible attempt to overcome this problem). 

We have seen in section 111 CI that the isocurvature 
mode contribution can create distortions in the small- 
scale linear matter power spectrum. Of course, this extra 
freedom was not taken into account in the definition of 
the grid of models in VHS. In principle, we should run 
simulations for a new grid with extra parameters (a, /3, 
"iso, "cor)- Alternatively, we can carry a tentative analy- 
sis with the same function bp(k) and the same error bars 
as in the pure adiabatic case, and check the validity of 
our results a posteriori. The idea is simply to select a 
marginally excluded model with the largest possible de- 
viation from adiabaticity in the matter power spectrum. 
For this model, we run a new hydrodynamical simula- 
tion and we compare P F {k)/P{k) with the function b 2 F (k) 
used in the analysis. In case of good agreement, the re- 
sults will be validated. We expect this agreement to be 
fairly good on large scales, but deviations should appear 
on small scales, because of the different non-linear evolu- 
tion. 

The use of state-of-the-art hydrodynamical simulations 
is a significant improvement compared to previous stud- 
ies which used numerical simulation of dark matter only 
H3. We use the parallel TreeSPH code GADGET-II [i| 
in its TrccPM mode which speeds up the calculation of 
long-range gravitational forces considerably. The simu- 
lations are performed with periodic boundary conditions 
with an equal number of dark matter and gas particles. 
Radiative cooling and heating processes are followed us- 
ing an implementation similar to |50l | for a primordial mix 
of hydrogen and helium. The UV background is given by 
|59j . To maximise the speed of the simulation a simplified 
criterion of star formation has been applied: all the gas at 
overdensitics larger than 1000 times the mean overden- 
sity is turned into stars [34| . The simulations were run on 



COSMOS, a 152 Gb shared memory Altix 3700 with 152 
CPUs hosted at the Department of Applied Mathematics 
and Theoretical Physics (Cambridge). 



C. Systematics Errors 

There is a number of systematic uncertainties and sta- 
tistical errors which affect the inferred power spectrum 
and an extensive discussion can be found in (33, H3, H3, 
VHS estimated the uncertainty of the overall rms 
fluctuation amplitude of matter fluctuation to be 14.5 % 
with a wide range of different factors contributing. 

We present here a brief summary. The effective opti- 
cal depth, r ff = — ln(F) which is essential for the cali- 
bration procedure has to be determined separately from 
the absorption spectra. As discussed in VHS, there is 
a considerable spread in the measurement of the effec- 
tive optical depth in the literature. Determinations from 
low-resolution low S/N spectra give systematically higher 
values than high-resolution high S/N spectra. How- 
ever, there is little doubt that the lower values from 
high- resolution high S/N spectra are appropriate and 
the range suggested in VHS leads to a 8% uncertainty 
in the rms fluctuation amplitude of the matter density 
field (see Table 5 in VHS). Other uncertainties are the 
slope and normalization of the temperature-density rela- 
tion of the absorbing gas which is usually parametrised 
as T = T (1 + <5b) 7_1 . T and 7 together contribute up 
to 5% to the error of the inferred fluctuation amplitude. 
VHS further estimated that uncertainties due to the C02 
method (due to fitting the observed flux power spectrum 
with a bias function which is extracted at a slightly differ- 
ent redshift than the observations) contribute about 5%. 
They further assigned a 5 % uncertainty to the somewhat 
uncertain effect of galactic winds and finally an 8% uncer- 
tainty due the numerical simulations (codes used by dif- 
ferent groups give somewhat different results). Summed 
in quadrature, all these errors led to the estimate of the 
overall uncertainty of 14.5% in the rms fluctuation am- 
plitude of the matter density field. 

For our analysis we use the inferred DM power spec- 
trum in the range 0.003 s/km < k < 0.03 s/km as given 
in Table 4 of VHS. (Note that, as in Q we have re- 
duced the power spectrum values by 7% to mimick a 
temperature-density relation with 7 = 1.3, the middle of 
the plausible range for 7 |6(|). 

Unfortunately at smaller scales the systematic errors 
become prohibitively large mainly due to the large con- 
tribution of metal absorption lines to the flux power spec- 
trum (see Fig. 3 of Ref. j34j) and due to the much larger 
sensitivity of the flux power spectrum to the thermal 
state of the gas at these scales. 
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07 8 m 1 o 10 ^ 20 55 65 75 FIG. 3: Two-dimensional likelihood for the amplitude of the 

isocurvature mode and of the cross-correlation component, 
near the pivot scale. We adopted a flat prior within the ellipse 

FIG. 2: Likelihood for the AD+CDI model, using all our ( which ap p ears here as a circle) in which these parameters are 

data set. The first eleven parameters are independent, while defined 

the last four are related parameters (with non-flat priors). 



parameter 


la C.L. 




0.0235 ± 0.0011 




0.125 ± 0.005 


e 


1.045 ± 0.008 


T 


0.11 ± 0.05 


n a d 


0.97 ± 0.02 




1.9 ± 0.5 


(5cor 


within prior range 


log[10 10 (yi 2 + B 2 )] 


3.3 ± 0.2 


a 


< 0.20 


2/3[a(l -«)] 1/2 


0.1 ± 0.2 


-^Ly — a 


0.8 ± 0.2 


ru 


0.68 ± 0.03 




0.88 ± 0.06 




13 ± 4 


H 


69 ± 3 



TABLE I: la confidence limits for the AD+CDI model, using 
all our data set, for the eleven basis parameters with flat 
priors, and below, for related parameters. 



IV. FITTING THE DATA 

A. Parameter basis and priors 

Any AD+CDI model is described by the usual six pa- 
rameters of the ACDM model, plus four parameters for 
the isocurvature sector (two amplitudes and two tilts). 
Like in most of the literature, we define the amplitudes 
parameters at the pivot scale ko = 0.05 Mpc -1 . For the 



isocurvature fraction, we could decide to impose a flat 
prior on /; so , or a, or any function of them; different 
choices are not equivalent, in general. We will come back 
to the dependence of the final result on the choice of pri- 
ors in section II V El Meanwhile, we chose a specific set 
of parameters which appear linearly in the expression of 
the observable power spectra, a and 2/3 -y/ a(l — a), and 
that we believe arc physically relevant. As already men- 
tioned, these two parameter are defined within an ellipse, 
in which we assume a flat prior. Furthermore, we must 
take into account the inequality 

|cosA| = |/3|(^) C "<1 (13) 

which should hold at least over the scales probed by the 
data, i.e. typically between fc m ; n = 4x 10 -5 Mpc -1 and 
&max = 2 Mpc -1 . This is achieved by introducing a new 
parameter <5 cor = n COI / In |/3| _1 , with a flat prior within 
the range —0.14 < <5 cor < 0.27. In summary, our basis 
parameters with flat priors consists of: 

• the baryon density, uih = f^b^ 2 , 

• the cold dark matter density, uj c = Q, c h 2 , 

• the ratio 9 of the sound horizon to the angular di- 
ameter distance multiplied by 100, 

• the optical depth to rcionization, r, 

• the adiabatic tilt, n a d, 

• the isocurvature tilt, ni so , 
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FIG. 4: Likelihood of the isocurvature-related parameters, 
for the three combinations of data sets described in section 
11V CI "Lyman-a" (red), "2dF bias prior" (green) and "none" 
(blue). (Left) Marginalized la and 2a confidence levels in the 
(a, 2/3[a(l — a)] 1 ^ 2 ) space. (Right) Marginalized probability 
distribution for ni so . 



• the parameter related to the tilt of the cross- 
correlation angle, 8 COI G [—0.14,0.27], 

• the overall normalization, ln[10 10 (A 2 + B 2 )], 

• the isocurvature fraction, a, 



• the cross-correlation amplitude, 2/3-y/ a(l — a). 

In addition, there are three independent parameters re- 
lated to observations: the Lyman-a calibration parame- 
ter ylLy- Q defined in |6lj . on which we impose the same 
Gaussian prior A^ y ^ a = 1.0 ± 0.29; and the two bias 
parameters associated to the 2dF and SDSS data with 
flat priors. Our full parameter space is therefore 13- 
dimcnsional. 



B. Results 

We compute the marginalized Bayesian likelihood of 
each parameter with a Monte Carlo Markov Chain 
method, using the public code CosmoMC . The re- 
sults are displayed in Fig. [2] and Table [i] (after marginal- 
ization over the 2dF and SDSS bias parameters). The 
data favors purely adiabatic models, but remains com- 
patible with an isocurvature fraction a < 0.40 at the 
2a (95%) confidence level (CL), with a tilt 7ij so = 
1.9 ± 1.0 (2<T CL). The one-dimensional likelihoods for 
a, 2/3 y a(l — a) must be interpreted with care: the 
fact that these parameters are defined within an el- 
lipse implies that there is more parameter space avail- 
able near a = 0.5 and 2/3y / a(l — a) = 0. More in- 
teresting arc the two-dimensional likelihood contours for 
(a, 2(3y/ a(l — a)) displayed in Fig. since in this rep- 
resentation the prior is really flat inside the ellipse. From 
this figure, it is clear that the data prefers an uncorre- 
cted isocurvature contribution. The flatness of the 5 cor 
likelihood shows that the data give no indication on the 
tilt of the cross-correlation angle. 
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FIG. 5: Favored ranges for the matter power spectrum 
P(k) in the three runs "Lyman-a" (dark), "2dF bias prior" 
(medium) and "none" (light), compared with our Lyman-a 
data, from the LUQAS quasar spectra (left) and from the re- 
analyzed Croft et al. spectra (right). The bands represent 
the envelope of all the matter power spectra in the Markov 
chains (after eliminating models with the worse likelihood). 
Each power spectrum has been computed at the median red- 
shift of the data and re-expressed in units of km/s. In addition 
to the statistical errors, the data points share an overall effec- 
tive calibration error, whose standard deviation is displayed 
in the top right corners. For the run including the Lyman-a 
data, each power spectrum has been divided by the value of 
the calibration parameter. The red dashed curves show the 
particular power spectrum discussed in section HV Dl 



C. Specific impact of the Ly-a data 

The Lyman-a forest provides a powerful indication on 
both the amplitude and the shape of the matter power 
spectrum for k > 0.01 s/km, i.e. roughly larger than 
1/i/Mpc. In order to illustrate the importance of this 
data set in our results, we repeat the same analysis with- 
out Lyman-a data. In this case, there are two options: 
we can either use the 2dF and SDSS galaxy power spec- 
trum data as a constraint only on the shape of the matter 
power spectrum, as already done in the previous analysis 
of section liV Bl or introduce a bias prior derived e.g. from 
the third and fourth-order galaxy correlation function of 
the 2dF catalogue [3^, , in order to keep an informa- 
tion on the amplitude of the matter power spectrum|67| . 

For these three cases, that we call "Lyman-a", "2dF- 
bias prior" and "none", the 2a upper bound on a are 
respectively equal to 0.4, 0.5 and 0.5. The likelihoods for 
the most interesting parameters are displayed in Fig. 0] 
As expected, the Lyman-a data set is significanty more 
powerful than the 2dF bias prior for cutting out mod- 
els with large a, and even more clearly, with large n\ sa or 
large anticorrclation, as can be seen in Fig. 4. It is impor- 
tant to note that without these data, all results depend 
on our arbitrary prior n\ SQ < 4: values far beyond this 
upper bound could still be compatible with the data, as 
also found in Ref. when using the same pivot scale. 
In the presence of the Lyman-a data, we get a robust up- 
per bound on nj SO , and none of our priors play a role in 
the final results, with the exception of the well-motivated 
<5 cor prior. 
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FIG. 6: The ratio b F (k) 2 = P flux (/fc)/P m atter(fe) at 2 = 2.75, 
computed from the hydrodynamical simulations as the ratio 
of the flux power spectrum (averaged over 1000 line-of-sights) 
over the input linear matter power spectrum. The solid blue 
curve shows the result for an adiabatic ACDM close to the 
best-fit model, while the dashed red curve was obtained from 
the "most extreme AD+CDI mixed model" defined in sec- 
tion [^^0] and here labelled as ACDI. The green band shows 
the region in which the Lyman-a data is used in the present 
analysis. 

The impact of the Lyman-a data can be understood 
visually from Fig. After running each case, we con- 
sider the collection of all matter power spectra in our 
Markhov chains (except models with a bad posterior like- 
lihood C < £ m ax/5). The gray bands in Fig. |S] cor- 
respond to the envelope of all these P(k)'s, compared 
to the Lyman-a data points. As expected, when the 
Lyman-a is not used, the band gets very wide above the 
wavenumber k ~ 0.2 h/Mpc ~ 2 x 10~ 3 s/km (note 
that for models with ni so = 4, the small-scale power 
spectrum is asymptotically flat). The role of the bias 
prior is marginal: it simply favors models with the lowest 
global normalization, but without affecting the isocur- 
vaturc fraction and tilt. Using the Lyman-a data, we 
can exclude any break in the power spectrum on scales 
k < 3 /i/Mpc ~ 3x 10~ 2 s/km. This results in much 
stronger constraints for the parameters (a, (3, ni so ), as can 
be seen from Fig. 4. 

D. Checking the validity of the Ly-a data fitting 
procedure 

We apply the strategy described in section ITlI Bl in or- 
der to check the validity of our Lyman-a data fitting 




0.5 10 0.5 1 



FIG. 7: Two-dimensional 68% and 95% confidence limits for 
the CDI mode amplitude and cross-correlation angle (eval- 
uated at the pivot scale). The two plots show the results 
of two independent runs with different parameter basis and 
priors. On the left (solid red curves), the parameters are 
(a, 2/3-y/ a(l — a)), with a flat prior within the ellipse. On 
the right (solid black lines), the parameters are (/i so ,cosA), 
related to the previous parameters through eqs.JHJ, with a flat 
prior within the square. The dashed curves show for compar- 
ison the likelihood contours obtained for one parameter set, 
assuming a flat prior on the other parameter set. 



procedure. We take the large number of samples con- 
tained in our Markov chains, and eliminate all models 
with a likelihood smaller than £ max /10 (in terms of ef- 
fective x 2 , this corresponds to A% 2 = x 2 — Xmin > 20). 
We then select the model with the largest value of a, 
which represents the strongest deviation from the purely 
adiabatic model. The corresponding matter power spec- 
trum is plotted in Fig. [S] and has a break around k ~ 
5 h/Mpc ~ 5 x 10 _2 s/km. Above this wavenumber, the 
slope of the power spectrum is given by eq. I|1(J|) with 
niso = 2.7. For this "extreme" model, we perform a 
hydrodynamical simulation as described in section Ull Bl 
and compare the bias function bp{k) with that assumed 
throughout the analysis. As shown in FigEI in the range 
0.003 < k < 0.3km/s probed by the data, the difference 
between the two functions is very small with respect to 
the statistical errors on the data. We conclude that in 
the present context, our Lyman-a data fitting procedure 
is robust, and does not introduce an error in the la or 
2cr bounds derived for each parameter of the AD+CDI 
mixed model. 



E. The role of parametrization and priors 

The fact of choosing a top-hat prior in the 
(a, 2/J-\/a(l — a)) parameter space is rather arbitrary. 
Other groups prefer to take top- hat priors on /j so , de- 
fined in Eq.©, and cos A = (3. Due to the non-linear 
transformation between the two basis, they are clearly 
not equivalent in terms of priors (see the discussion of 
this point in [3l]], in the context of Baycsian Evidence 
calculation for adiabatic versus mixed models). 

We checked this issue explicitly with an independent 
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FIG. 8: Likelihood for /i SO in the curvaton model. 



run based on the (/i SO , cos A) basis. The results are sum- 
marized in Fig. As expected from the Jacobian, the 
(/iso,cosA) option gives more weight to models with a 
small isocurvature fraction. For instance, the run with a 
flat prior on /i so gives a la bound /; so < 0.26, while that 
with a flat prior on a gives /i so < 0.66. However, at the 
2a level, the relative difference is small (/i so < 0.75 ver- 
sus /iso < 0.87) because the Jacobian is asymptotically 
flat. 

In principle, in Baycsian analysis, the choice of param- 
eter basis and priors should reflect one's knowledge on the 
model before comparision with the data. However, in the 
absence of a unique underlying physical model motivating 
the presence of isocurvature modes, different scientists 
might put foward different choices of prior. This intrin- 
sic freedom in Bayesian analyses should always be kept 
in mind when quoting bounds, especially for parameters 
which represent physical ingredients not strictly needed 
by the data, which is the case here for the isocurvature 
sector parameters (for other parameters such that the 
data picks up a narrow allowed region, a change of pri- 
ors won't affect the bounds very much). However, even 
for the isocurvature parameters discussed here, it is re- 
assuring to see from our analysis that the 2a contours 
obtained from the the two runs and compared in Fig. 
are roughly in agreement. 



F. The curvaton model 

In this section we derive bounds on the specific case of 
curvaton models. The curvaton hypothesis is an ingen- 
uous way to generate the observed curvature perturba- 
tion from a field (the curvaton) different from that which 
drives inflation (the inflaton) 0|. In practice there is 
not much difference in the phcnomenological signatures 
left in the CMB and LSS compared to an ordinary infla- 
tionary model. However, there are a few cases in which 
it is possible to leave a "residual" isocurvature compo- 
nent, together with the dominant curvature contribution. 
More specifically, in the curvaton models in which the 
curvaton field is responsible for the CDM component of 
matter, there are various possibilities depending on the 
time of creation of CDM versus the decay of the curvaton 



field. In all these cases, the curvature and isocurvature 
perturbations are related to the gauge invariant Bardeen 
variable £ as 

S = 3(Ccdm - C) , (14) 

K = -(. (15) 

Let us classify here the different cases: 1) when CDM- 
creation occurs before the curvaton decays and the frac- 
tion r of the total energy density in the curvaton field at 
the time of its decay is negligible. Then S = — 3C = 37tL, 
which corresponds to /j SO = 3 (a = 0.9), and (3 = +1 
(maximally correlated), with n; so = n a( j; 2) when CDM- 
creation occurs before the curvaton decays but the frac- 
tion r at decay is important. This case requires specific 
model input and in principle can have any value of /i SO 
and f3, while nj SO = n a d; 3) when CDM-creation occurs at 
the decay of the curvaton and the fraction r < 1. In this 
case, Ccdm = C/ r an d tnus <-> = 3(r — 1 — 1)C = 3(1 — l/r)1Z, 
which corresponds to /; so = 3(1 — 1/?'), i.e. (3 = — 1 (max- 
imally anticorrclated) and rii so — n a( j; 4) when CDM- 
creation occurs after the curvaton decay. Then there 
is only one thermal fluid in equilibrium, Ccdm = C> an d 
there is no way to generate an isocurvature perturbation, 
5 = 0. 

Since case 1) is already excluded at many sigma, and 
case 2) is essentially identical (except for n; so = n a ^) to 
our generic analysis, we will concentrate on case 3) of 
a maximally anticorrclated mixture of isocurvature and 
adiabatic modes with equal tilts and 5 cm = 0. Our results 
are summarized in Fig. [8] which shows the likelihood dis- 
tribution for the generic curvaton model. We have used 
"ad = "iso, ^cor = 0, and f3 = ±1, which is equivalent to 
f3 = 1 and /; so positive or negative: /i so > corresponds 
to /3 = 1, or positive correlation between lZ ra d and 5 ra d, 
i.e. suppression of power in P(k) and in the large-scale 
CMB temperature spectrum; while /i SO < corresponds 
to the opposite anti-correlated case. 

We find / iso = 0.04±0.09 at the 2cr-level, which implies 
r > 0.98 at the same CL. In our opinion, such a stringent 
constraint on the fraction of energy density in the cur- 
vaton at decay calls for a tremendous finctuning (there 
is no physical reason to expect that the curvaton should 
decay precisely when it is starting to dominate the total 
energy density of the universe, within 2%), which makes 
the curvaton hypothesis in its most attractive scenario 
very unlikely. 

G. The double inflation model 

Another chance to generate an observable isocurvature 
signature is through the pos sible presence of two scalar 
fields driving inflation 0, l65j . The simplest case at hand 
is that of two massive fields coupled only gravitationally: 

£ = \{<t>hw<t>h - m l<t>h) + -^i-^T - > ( 16 ) 
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where ttj/j and mi are the masses of the heavy and light 
fields respectively. 

We assume slow-roll conditions during inflation, and 
use the number of e-folds till the end of inflation s = 
— ln(a/a en d) to parametrize the fields as: 



& = 2^G Shl2 6; ^ = 2^G C0S2 6 ' 

Using the field and Friedmann equations, we can solve 
for the rate of expansion during inflation: 



H 2 {s) ~ -s ■ mf [1 + (R 2 - 1) sir 



(18) 



where R = m^/mi, and find the number of e-folds as a 
function of 9: 



s(8)=s Q 



( sin g)2/(fl 2 -i) 
cos60 2fl2 /(« 2 -i) 



(19) 



The perturbed Einstein equations can be solved for 
long wavelength modes in the longitudinal gauge. As- 
suming that the heavy field decays into CDM whereas 
the ligth field produces other species, we find the magni- 
tudes of the curvature and entropy perturbation at hori- 
zon crossing. During radiation domination and for super- 
Hubblc modes, this gives: 



ftrad(fc) = - 
<Sra,d(k) = 



4ttG 



tj 1/2 



4ttG _i/2 

~W HkSk 



sin 9 k e h (k) + cos 9 k e ; (k) 
e fc (k) R 2 ei (k) 



sin t) k 



cos k 



(20) 



where e^(k) are gaussian random fields associated with 
the quantum fluctuations of the fields, and the subindcx k 
implies the value of the corresponding quantity at horizon 
crossing during inflation. One typically expects s k — 60. 
It can be seen from that the correlation power spec- 
trum has no scale dependence, and thus, for this model 
"cor = 0, while the adiabatic and isocurvature tilts have 
expressions 

{R 2 - 1) tan 2 6 



"-ad = 1 

Sk 



2s k (l + R 2 tan 2 9) 2 



(21) 



(R 2 - 1){R 2 tan 4 6 - 1)(1 + tan 2 9) 
Ulso ~ 1_ Sk (l + R 2 t<m 2 9) 2 (l + RHtm 2 8) ' (22) 

whose values, for s k = 60, are typically rt a <j = 0.97 and 
n iso in the range [0.97,0.90] for R e [1,4]. Since 7i iso > 
0.93 at 95% c.L, models with large values of R are ruled 
out. 

It was shown in |2l| that a relationship beteween a and 
(3 can be found. It can be simply expressed as a straight 
line in our parameter space: 

2/Va(l-<»)= 2(fl2-1) (l-«)- (23) 

Sk 



On the other hand, for these models, the parameters a 
and P have minimum and maximum values respectively, 
which only depend on the ratio R and the number of 
c- folds Sk, 



On 



(R 2 



I) 



4 + (R 2 

R 2 - 1 



2/V"(l-a) 



R 2 + l ' 
2s k (R 2 



I) 2 



1) 



+ (R 2 + l) 2 



(24) 
(25) 
(26) 



Using the results of section lTV Bl we find that the inclu- 
sion of the Lyman-a data significantly improves the pre- 
vious bound on R to R < 3 at 95% c.l. This bou nd comes 
mainly from a combination of bounds on 2/3-^/ a(l — a) 
and n iso . 

We did not find necessary to generate a n COI = sam- 
pling for this model. In our results, the parameter S COT 
has a flat distribution and thus is unconstrained. We 
therefore expect similar results when fixing it to zero. 



V. CONCLUSIONS 

In addition to CMB, LSS and SNIa data, we used 
some recent Lyman-a forest data to further constrain the 
bounds on possible CDM-isocurvature primordial fluctu- 
ations. We find that the systematics induced — in partic- 
ular, those associated with the recovery of the linear dark 
matter power spectrum from the flux power spectrum — 
are greatly compensated by the valuable information on 
the small-scale matter power spectrum provided by the 
Lyman-a data. 

Before summarizing our results, it is worth mentioning 
that when we omit the Lyman-a forest data our bounds 
agree very well with those of Ref . ■ The authors of 
work with a pivot scale ko = 0.02 Mpc" 1 , but they also 
show how their results arc modified when they take ko = 
0.05 Mpc -1 like in the present paper: in that case the 
agreement with us is particularly good. The comparison 
of our results with the WMAP analysis from Ref. |24| is 
more puzzling: using or not some Lyman-a data, they 
always find much stronger bounds on /; so than us. It is 
true that we have one more free parameter S COI , and that 
we do not introduce a prior on the 2dF bias; however, 
even when we fix <5 cor = and introduce such a prior, 
our /i S o bound remains much weaker. So far, private 
communications with the WMAP team did not allow us 
to understand the origin of the discrepancy. 

Using all our data set, we find at the 95% confidence 
level, an isocur vature fra ction a < 0.4, a cross-correlation 
amplitude 2/3 ^/a(l — a) = 0.1 ±0.4, and an isocurvature 
tilt Tii so = 1.9 ± 1.0. The tilt of the correlation angle 
remains unconstrained. If we switch to the basis used for 
instance in Ref. [24| we find /i so < 0.75 at 95% c.l. 
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In the case of a curvaton scenario where CDM-creation 
occurs at the decay of the curvaton - a case in which the 
adiabatic and isocurvature modes are maximally anti- 
correlated, (3 = — 1, and n a d = n- lao - we find /; so < 0.05, 
still at the 95% confidence level. This requires that the 
fraction r of the total density in the curvaton field at that 
time be fine-tuned between 0.98 and one. Finally, if we 
assume a double-inflation model with two massive infla- 
tons coupled only gravitationally, such that the heaviest 
field decays into CDM, while the lightest one into stan- 
dard model particles, we find that the mass ratio should 
obey R < 3 (95% c.L). 
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